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SIWMARI 


The problem, of "the ■bhree-polnt body of revolution that has Ttilnlmum 
■vreive drag; baaed, on linearized, stper sonic-flow theory, has been solved, 
formally withotcb resort to slender-body approximations. The evaluation 
of the source-distribiation function, els well as the body shape, requires 
numerical procedLures. Conplete calculations are made for the special 
case of a closed body with, a given section (the third point) at the mid- 
body position. The existence of a comer at the third point is found to 
be characteristic of a li^ar-theory solution of the minimum-drag three- 
point problem. The conclusion is reached that the influence between the 
forward anil reaiTward portions of a closed body results in a sitial.l increase 
in the radii of the minimum-drag shape. 


BITRODUCnON 


A large number of minimum-drag problems have been solved by using 
the slender-body approximation to the .linearized supersonic-flow theory. 
In 1935 Von KaTTnan (ref. 1) determined the projectile tip that has mini- 
mum wave drag according to slender-body theory. Later Sears (ref. 2) 

ATifl Haack (ref. 3 )# using slender-body theory, determined minimum-wave- 
drag shapes for projectile tips and closed bodies of revolution subject 
to various combinations of auxiliary co n d it ions of constant length, con- 
stant caliber, and constant volume. 

Ferrari (refs. 4 and 5 ) has considered the minimum-drag problem for 
the length- caliber body and the ducted body on the basis of linear theory 
without resorting to the slender-body approximation. Because he used 
the linearized drag integral in the triple-integral form, Ferrari found 
it necessaiy to \ise a series expansion for the source-distribution func- 
tion, "sdilch introduced considerable con^lexlty at an early stage in the 
analysis. Recently Parker (ref. 6) reduced the linearized drag integrea 
for bodies of revolution to a dovible-integral- form and determined, with- 
out resort to slender-body theory, the minimum-wave-drag shape for a 
transition section between two semi-infinite c yl i n d er s, a special case 
of which is the projectile tip. 
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This paper presents a linear-theory solution to the problem of deter- 
mining the three-point body of revolution that has minimum wave drag. Q?he 
body configuration is the transition section of reference 6 with the a d di- 
tional requirement that the section must have a given arbitrary radius at 
a given arbitrary position (the third, point) . The linearized, drag integral 
in double-integral form from .reference 6 is used, and two auxiliary condi- 
tions similar to the one utilized in reference 6 are employed. The soTxrce 
distribution for mimmmti drag is obtained in a closed form involving a 
single integreil of elliptic type idiich, presumably, is most conveniently 
treated numerically. Determination of the shape requires an additional 
numerical procedure. A closed body with the maximum section at the mid- 
point is chosen for a conplete calculation. The minimum-drag shape has 
a comer at the position of the third point. For this special case the 
shape is slightly asymmetric. The radii downstream of the maximum sec- 
tion are abo\xt 1 percent larger than the corresponding reidii upstream. 


SYMBOIS 


a axial coordinate of fixed third point 

C constant factor in source-distribution function for special 

body, ft /sec 

D drag, lb 

f source-distribution function 

f ' derivative of f with respect to its argument 

L an axial distance in equation (l) such that when x ^ L there 

is no momentum flow through the control cylinder of radius 

M free-stream Mach number 

r radial distance in cylindrical coordinates 

R radius of general point on body 

Rq radius of upstream cylinder 

R^ radius of downstream cylinder 

Rg. radius of fixed third point 

_2 step function (\mity for negative argument, zero for positive 

argument) 
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U free-stream velocity, fb/sec 

Vp radial velocity, ft/sec 

X axial distance in cylindrical coordinates 

p = \/m2 - 1 

5f ' variation in derivative of source-distribution function 

p free-stream density, slxigs/cu ft ‘ 

7vi,?V 2 Lagrangian multipliers 

|,Ti dummy variables of integration 

a axial cooidinate -where sovcrces begin 

All distances are made dimensionless by using the length, of the body 
as the unit of measurement. 


ANAI2SIS 

The linearized drag Integral for a body of revolution is (ref. 6) 


D = jtp 


y^L-pRl j 


L-PRl 

f'(l) f‘(n) cosh ^ 


(L - |)(L - Ti) - 


^Rl(l - k) 


d| dTi 


( 1 ) 

Equation (l) gives the drag of an axial distribution of sources f(|) 
which begins at a and is svibject to the condition that there be no 
moment-um flow through the control cylinder of radius for x > L. 


Figure 1 is a schematic representation of the general configuration. 
The do-wnstzeam and upstream cylinders are of radius R^ and Rq^ respec- 
tively. The body contour is required to have -the radius R^ at x = a, 
the third point of the "three-point problem." The distances Ra and a 
are gi-ven, arbitrary valvies; that is, the drag is not minimized -with 
respect to the location of the third point. Thd radius of the control 
cylinder of equation (l) is tahen els Ri, the radius of the downstream 
cylinder. For all lengths, the length of the body is xised as the unit 
of measinrement . 
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The two aTJxlIlary conditions require that the point (l^Rl) be 
on the same streamline as (0,Bo), that is. 


pTSt(Ri2 - Bq 2) = 2i£p 


r 

d oft 


^(El-Ro) 


HlVr(x,Ri) dx 


( 2 ) 


and that the third point (a,Ra) the same streamline as (0 ,Rq), 

that is. 


pa 


pIit(Ra^ - Rq^) = 2itp / RaVr(x,Ra) dx 

-'3(Ra-Ro) 


(5) 


Erom the general solution of the linear st^ersonic-flow equation for 
the case of axial symmetry, the radial velocity is given by (ref. 1) 


Vr(x,r) 


i f»(S)(x-£)d| 

^ J.« ^(x - 1)2 - p^2 


(4) 


In the present problem f’(|) vanishes for | < -^Ro = a. Inserting 

equation (4) into eqoxations (2) and (3), interchanging the order of inte- 
gration, and integrating over x yiel^ 


2(Ri2 - Rq 2) = f ^ f • (I) \/(rr7)2T^^^ d| (5) 

2 ^ -PRO V 


and 



- Rq 2) = r f'(l)\/(a - g)2 - p 2 Bg ^2 ^ 

d _pEo V 


(6) 
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For the minimizing process -vdiich is to follow, it is mathematically 
convenient to change the Integral in equation (6) , to an Integral over 

the range -pEo to 1 - pR^ by defining a step function _Z(l-a+PRa:) 
such that 


l(&-a+pEa) = 1 (I ^ a - pRa) 


^(l-a+pRa) = 0 (I > a - pRa) 


Therefore, equation (6) may be written 


n 1-pRn I 

g(Ra2 - V) = / i(l-arfpRa) f(l)\/(a - |)2 - p^^^ d| (8) 

2 J V 


The problem of minimizing the drag given by eqiiation (l) with 
a = -PRq and L = 1, subject to the auxiliary conditions given by equa- 
tions (3) and (8), can be esspressed in the usual manner of the calculus 
of variations. The e3q>ression 

o f(S) f(n) co8h-i (J- - - pV at ai, + 

J-PBo PtLd - n) 


f'(t) iflTt)2T^^at + 2?i2 r ^J(t-wpBa) f'{t)\/(a - t)^ - 
J-PRO '' '^-PBO 


(■vdiere constant factors have been absorbed into the Lagranglan multipliers 
?vi and A2) must have a zero variation corresponding to an arhitrary RTtial l 

variation 6f* in f*; that is, 

G(f’+6f',a,Ro,Ri,Ra) - G(f ’,a,Ro,Ri,Ra) = 0 (lO) 


Combining two terms, after interchanging dummy variables in one, and 
neglecting the term involving (8f')^ yields 
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‘l-pRi 


f ' " 8f‘(|) 

^-PBO 


p 1-PR]_ 

/ f'(Tl) cosh."^ 

-'-PEO 


(1 - I) (1 - Ti) - p2Ri2 


pRiCl - Tl) 


dT) + 


7^1 ^(1 - 1)2 - p%i2 + 7^2 i(S-a+pEa) ^(a - |)2 - p^^ 


d| = 0 


( 11 ) 


Since the variation 5f'(S) is arbitrary. 


,1-pRi 


f(.i) cosh-1 

-'-PEO 


(1 - |)(1 - il) - P%i2 


pRiCl - n) 


dT] + ^i\Ju - 6)^ - P^l^ + 


?^ 2 i(l-a+pBa)^ (a - l)^ - p^r/ = o 
An integration by parts, with f (-pEo) = 0, yields 


(12) 


/ 


1-pEi 


^(ti) dTi 


^ j/(l - tj) 2 . p%^2(| _ 


1/(1 - 1)^ - p%^2 


= H(l) 


(15) 


or, if 


then 


f(^) 


^(1 - ti )2 - p 2 R ^2 


= P(^) 


/ 


-PEo I - n 


(14) 
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The solution of equation (l4) is 


F(l) = 


K - 


,t2j^(l - pRi - |)(| + PEo)L 

\/(l - pRi - Tl) (n 4- PEq) H(ti) dTi 


-PKq 


I - n 


(15) 


where K is a constant and' | is in the range -PEq =6 = 1- ^Rl. A 

discussion of the conditions for th.e existence of this solution is given 
by Tricomi in reference 7* 

It is a straightforvraurd procedure to show that the solution of equa- 
tion ( 13 ) is 


f(l) 




PBl - IJ 


^B0 


|K + 


(- 


- - PEq 




oa-pRa 

(a + pRa - n) (a “ PRa - n) (^ + pRa) 

1/2 

dTl 

^-PRO 

1 + pRi - T| 

T - 6 


( 16 ) 


A convenient form of the solution, which can be obtained by algebraic 
manipulation, is 


f(l) 


= 1 A 


pBi - I 


A 2 


PRo 

l(l4fRi) - 1(g) 
1 + pRl - I 


K + JtAqI 




- P^i - Pfip 


->) 


( 17 ) 


where 


1(6) = r 

- 


a-pRa 

PRO 


(a + pRei - Tj) (a - pRa - ti) (t) + PRq) (l + PBi - 



dTi 


I - n 


(18) 
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Requiring that f(-pEo) = 0 yields 


+ pRi - 6)(l + PBo) + >2 


f(l) 


1 + pRl - E 


1 + pH.j_ + pEo 


ji(-pHo> - l(l+pRi)] + l(l+pRi) - 1(6) 


it2^(l + pEi - 6)(6 + PBo) . 


(3S) 


The constants Aq and A 2 are determined by satisfying the auxiliary 

conditions, equations ( 5 ) and ( 6 ) . The auxiliary conditions are special 
cases of the general shape expression for the points (a,Ro) and (l^Ri), 
and in the general case numerical procedures are necessary in order to 
satisfy them because l(|) could not be integrated. If in equation ( 19 ) 
Ag is set equal to zero, the result is the minimum-drag source- 
distribution function for the “two-point” problem of reference 6 . 

The general shape expression is 


■ (20) 
J-PKO . j)2 . p%2 

(which is eq. ( 27 ) of ref. 6, after interchanging the order of integra- 
tion and performing one integration) . The shape satisfies the Gothert 
similarity rule (ref. 8) so that one calculation gives a functional rela- 
tion between pR and x for given values of a, PRq, pR^^, and pRa 

and therefore gives a family of minimum-drag shapes. 


CALCULATIOMS FOR A CILOSED BOUT 


For a closed body, pRo = pR^ = 0, the source-distribution function 
(eq. ( 19 )) reduces to 


f(l) = 


it2^|(l - I) 


Jt7^i(|2 - I) + A2l[l(l) - 1(0)] + A2 [i(0) - I(0] 


( 21 ) 
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where 


pa-pRa r 1/2 

lU) = / (a + pRa - n)(a - pRa - - tj) - 23 — (22) 

The avixillary-condl'tion equation ( 5 )# after an integration by parts anfl 
use of f(o) = 0 , reduces to 



which, when use is made of tlie fact that 

° \/l(l - !) 

reduces the source-distribution function, equation (21), to 


f(|) = c - ^&(Q) + 1(1)] ( 23 ) 

\/l(l - I) 


where 


C = 



Since the last avudliary condition to be satisfied (eq. (6)), which, 
fixes the constairt 7^2^ is a special case of the shape e3q>ression, it is 

convenient to rewrite equation (20) for the closed body as 


-L(pr)2= f(s) (x - t)d{ ( 2 ^, 

Jo 0 . 5)2 . J8E2 


The value of f(|)/c 
of equation (24) for 


may be calculated from equation ( 23 ) • Evaluation 
X = a and pR .= pRa ( corresponding to the last 
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axtsdLliary condition) fixes the value of u/23^. Then equation (24) is 
a functio na l relation between pR and x giving the shape of the minimum- 
drag closed body. 

The complete calciilations were made for the case of a = 0*5 'and 
pRa = 0.2, a closed body with the third point at the symmetrical midpoint 

position. Details of the calculation are given in the appendix. 


DKCtBSION OP REBUEflS 


The problem of the three-point body of revolution that has minimum 
wave drag, based on linearized s\ipersonic-flow theory, has been solved 
formally. Equations (19) and (20) permit the calculation of the shape 
of the general three-point body that has miniTnum wave drag. Except for 
a slightly increased complexity in the first two integrations of equa- 
tion ( 20 ), corresponding to the two auxiliary conditions, the ceilcula- 
tions are essentially the same for a ducted body as for a closed body. 

An actual ducted-body shape was not computed. 

Figure 2 shows the integral l(|) and the quantity f(|)/c for 
the computed closed body (a = 0.5 and pRa = 0.2). The one-half -order 
singularity in the slope of l(|) and f(|) at | = 0.5 corresponds 
to the comer at x = 0 . 5 . The shape of the body is given in figure 5 
and table I. An inspection of the mathematics shows that the comer is 
present for any finite value of pRa* In the slender-body limit (pRa 

approaching zero) the comer vanishes as 2 (icpRa)^* 

Recent unpublished work by Clinton E. Brown of the Langley 
Aeronautical Laboratory shows that to the order of accuracy of linear 
theory the wave drag of bodies of revolution is reversible. With the 
plausible assxmption that the minimum-drag problem has a iinique solu- 
tion, the body compvited here should be symmetrical. The slight asym- 
metry fo\md (see table I) is not to be interpreted as a refutation of 
the reversibility theorem, since higher order terms, of appreciable 
Influence because the body chosen is rather thick, undoubtedly are handled 
(included or omitted) in different ways in the different treatments. 

Figure 4 permits an easy comparison of the minimum-drag body com- 
puted in this report with other minimum-drag bodies which have been com- 
puted. The ordinate is (Hbody " ®cone)/^Wix ®body radius 

of the body, Rcone radius of the cone whose vertex is at the nose 

of the body and whose base is the maximum cross section of the body, and 
Rpinv is the radivts at the TnayiimiTn section. The abscissa x is distance 

eilong the body and the maximum section is at x = 0.5* figure shows 
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(a) the forward and (B) the rearward portions of the body computed herein, 
(C) the coiresponding linear- theory minimum-drag projectile tip, (D) the 
Sears-Haack symmetrical length- caliber body, and (E) the Von TTp-rmpTi pro- 
jectile tip. 0?he radii of the minimum-drag bodies according to slender- 
body theory (d and E) are functions of x only, whereas the radii of the 
m i nimum- drag bodies according to linear theory (A, B, and C) are functions 
of pRa(=pikax) as well as x. 

It is apparent that the lack of symmetry of the body computed herein 
(the difference between curves A and B) is of the same order as the dif- 
ference between that body and two llnear-theoiy ml niTnTnn- fl-r q g projectile 
tips placed base to base. Both forward »nd rearward portions are thicker 
than the projectile tip and their average is thicker than the projectile 
tip In approximately the same proportion that the Sears-Mack length- 
caliber body is thicker than the Von Karman projectile tip. In the 
exbreme case where pEa. = O.5 (of course, coinpletely outside the range 

of validity of l ine ar theory) the body confuted in this report would 
become two Mach cones and the ctnrves A, B, and C would collapse into the 
abscissa. In the limiting case, pHa — ^0^ "the three-point body becomes 

the Sears-Haack length-caliber body; that is, c\nrves (A) and (b) coalesce 
Into curve (d) \dille curve (C) approaches curve (E) . 

Therefore, In the range of validity of linear theory, the minimum- 
drag length-caliber body with fixed caliber in the symmetrical position 
is thicker in both forward and rearward parts than the corresponding 
minimum-drag projectile tip. Thtis the concloision is reached that the 
Influence between the fonrard and rearward portions of a closed body 
results in a small Increase in the radii of the minimum-drag shape. 


CONCLUDING REMARKS 


The problem of the three-point body of revolution that baa Tnlnlnmin 
wave drag, based on linearized sv^ersonic-flow theory, has been solved 
formally without resort to slender-body approximations. The source- 
distribution function Involves an Integral which apparently cannot be 
evaluated in terns of simple functions and for which a series expansion 
may be used. Numerical quadratures are required to determine the shape 
of the body. 

The complete calculations are made for a closed body with a given 
section pRa =0.2 at the midbody position -a = 0.5. The shape is 
slightly asymmetric. A comer at the fixed-caliber position is found 
to be characteristic of a linear-theory solution of the minimum-drag 
length- caliber problem. 
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A coniparlsoii of "thfi confuted TnIniTmim body with other TtilTiltmnn bodies 
leads to the condusion that the influence between the forward and rear- 
ward portions of a closed body results in a small increase In the radii 
of the Tnj.nimum-drag shepe. 


Langley Aeronautical Laboratory^ 

National Advisory Committee for Aeronautics, 
La ngl ey Field, Va., March 15, I956. 
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APEEMDIX 


DETAILS OF THE CALCULATIONS 


The source of difficulty in maklug the calculations is the integral 
l(l) (eq.- (22)), which appears in the source-distribution function. 

For a = 0-5 aini pRa = 0.2, 


1 ( 6 ) = 



•7 - tl)(0.3 - n)Ti(l - 



(Al) 


C h anging the variables so that = 0.15(l + y) anfl p = 0.15(1 + t) 
yields 


[[0.15(1 + t)] = 0.15 \/(0.55) (0.85) 




fl-iy) 


U- } 

V 17 / 


-. 1/2 


\/l - y2 dy 
t - y 


(A2) 


If -l^t$l (O^l^ 0 . 3 ) the Integral is singular; if 1 < t ^ 17/3 

( 0.3 < 6 ^ 1) it is not singular. Since singular and nonsingular inte- 
grals of the type 



(n = 0, 1, 2, . . .) (A3) 


may be evaluated (singular integrals of this type may be deduced from 
Trlcomi’s paper, ref. 7)^ "ttie quantity 



was expanded in a power series in y. In this manner two series were 
obtained for l(|), one series for the range 0 ^ | ^ O .3 (-1 ^ ^ 1) 

and one series for the range 0.3 ^ 6 ^ 1*0 (l ^ t ^ 17 / 3 ) • Th.e com- 
putations were performed on a desk calculating machine and seven digits 
were carried. At most, twelve terms were required for seven-digit 
accuracy. Values of l(|) were confuted at approximately 120 points in 
the range 0 ^ | ^ 1, with the intervals smaller in the regions where 
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l(|) or f(g) changes rapidly. Equation (23) -was then used to calculate 
f(|)/c at the points at vhich l(|) had been evaluated. 

In the shape expression (eg.. (24)) 



The integral was evaluated either by Simpson's rule for graphical inte- 
gration or by replacing f(|)/c by a quadratic in | over a small range, 
depending iq)on how rgpidly the integrand was changing. First, equa- 
tion (a4) was evaluated for x = O .5 and ^R = 0.2, which fixed the 
value of the constant u/2p2c. For a typical point on the body, a value 
of X was chosen and the two members of equation (a4) were calculated 
for trial values of pR, and by Inteipolatlon a value of pR was foimd 
to satisfy equation (a4) . 

Figure 2 shows the variation of the integral l(|) and the quantity 
f(|)/C with Figure 3 and table I give the shape of the computed 
closed body, a = O.5 «rid pEg, = 0.2'. It is believed that the overall 
error of the numerical procedure is such that the fifth decimal place 
for pR is significant. 
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TABLE I 


CCaiPlICED ORDINATES FOR TEE CLOSED BODI WITH 


a = 0.5 AND 3Ea « 0.2 


X 

PR 

0 

0 

.1 

.05764 

.2 

.10120 

•3 

.13862 

.4 

.17141 

.49 

•19735 

•5 

.20000 

.51 

.19752 

.6 

.17283 

•7 

.i 4 o 47 

.8 

.10288 

.9 

.05866 

1.0 

0 







& 


Figure 2.- Tbe Integral l(|) anfl. the minlmum-cLrag source -distribution 
function f(6)/C. a = 0.5) P^a “ 0*2. 
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Oistonce along body, x 


Figure 4-.- Variation of ^4x)dy ” ^one)/^Wx 'wi'fc^ distance along "body 

for (a) the forward portion and (b) the rearward portion (x — >-(l - x)) 
of the body computed in this report, (c) the corresponding linear-theory 
projectile tip, (d) the Sears-Haach l en gth-callher hody, and (e) the 
Von Earidin projectile tip. 
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